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Using mainly numerical methods, we investigate the ground-state phase diagram of the S~2 
quantum spin chain described by H = J2 j i S j S J+i + S j S j+l+ AS I S i+i)+ D E^SJ) 2 , where A 
denotes the XXZ anisotropy parameter of the nearest-neighbor interactions and D the on-site 
anisotropy parameter. We restrict ourselves to the case of A > and D > for simplicity. Each 
of the phase boundary lines is determined by the level spectroscopy or the phenomenological 
renormalization analysis of numerical results of exact-diagonalization calculations. The result- 
ing phase diagram on the A-D plane consists of four phases; the XY1 phase, the Haldane/large- 
D phase, the intermediate-D phase and the Neel phase. The remarkable natures of the phase 
diagram are as follows: (1) there exists an intermediate-D phase which was predicted by Os- 
hikawa in 1992; (2) the Haldane state and the large-D state belong to the same phase; (3) the 
shape of the phase diagram on the A-D plane is different from that considered so far. We note 
that this is the first report on the observation of the intermediate-D phase. 

KEYWORDS: S = 2 quantum spin chain, ground-state phase diagram, quantum phase transition 



Quantum spin chain systems have been attracting in- 
creasing attention in recent years because they provide 
rich physics even when models are rather simple. The 
most well-known example may be the existence of the 
Haldane phase 1, 2 -* in integer spin chain systems only with 
isotropic nearest-neighbor (nn) interactions. 

In this study, using mainly numerical methods, we ex- 
plore the ground-state phase diagram of the S — 2 quan- 
tum spin chain described by the Hamiltonian 
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(1) 

a = x, y, z) is the a-component of the S — 2 
operator at the j-th site, and A and D are, respectively, 
the XXZ anisotropy parameter of the nn interactions 
and the on-site anisotropy parameter. Hereafter, we de- 
note the total number of spins in the system by N, as- 
sumed to be even, and the z component of the total spin 
by M(= £,•<?!)• 

The ground-state phase diagram on the A-D plane of 
the same Hamiltonian in the S=l case was discussed 
by Schulz 3 ) and den Nijs and Rommels, 4 -* and numeri- 
cally determined by Chen et al. 5 -* In the present S = 2 
case, Schulz 3 - 1 discussed the ground-state phase diagram 
with six phases obtained by the bosonization method; 
the phases are, in our terminology, the ferromagnetic 
(FM) phase, the Neel phase, the XY1 phase, the XY4 
phase, the Haldane phase and the large-D (LD) phase. 
The XY1 state is characterized by the power decay of 
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Fig. 1. Valence bond pictures of (a) the Haldane state, (b) the ID 
state and (c) the LD state. Large circles denote 5 = 2 spins and 
dots S =1/2 spins. Solid lines represent valence bonds (singlet 
pairs of two 5=1/2 spins, (l/yf^t-l - 4_t)). Two 5=1/2 spins 
in rectangles are in the (5tot , S* ot ) = (1,0) state and four 5 = 1/2 
spins in squares are in the (5tot, 5^ ot ) = (2, 0) state. 



the spin correlation function G±i(r) = (S^~ S 1+r ) and the 
exponential decay of G±4(r) = ({Si)} 4 {Si +r } 4 ), whereas 
the XY4 state by the exponential decay of G±i(r) and 
the power decay of G±4,(r). The valence bond pictures 
of the Haldane state and the LD state are shown in 
Figs. 1(a) and (c), respectively. About twenty years ago, 
Oshikawa 6 ) predicted, in .!5>2 integer quantum spin 
cases, the existence of the intermediate-D (ID) phase, 
the valence bond picture of which is depicted in Fig. 1(b). 
Figure 2(a) is an interpretation of Oshikawa's prediction 
by Aschauer and Schollwock 7 ) in the S — 2 case. 

On the other hand, carrying out the density-matrix 
renormalization-group (DMRG) calculation, Schollwock 
et al. 8 ' 9 ) and Aschauer and Schollwock 7 ) proposed the 
phase diagram in Fig. 2(b) and concluded the absence of 
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Fig. 2. (a) The shape of the ground-state phase diagram of the 
Hamiltonian (1), predicted by Oshikawa 6 ) and interpreted by As- 
chauer and Schollwock. 7 ) The symbols FM, XY1, XY4, LD, ID, 
H and N stand, respectively, for the ferromagnetic, XY1, XY4, 
large-D, intermediate-D, Haldane and Neel phases. The bound- 
ary lines between the XYl phase and the ID or H phase are 
not necessary to be A = 0. In other words, the XYl phase may 
intrude into the A > 1 region, (b) The shape of the ground-state 
phase diagram of the Hamiltonian (1), proposed by Schollwock 
et al. 7 ~ 9 > by the DMRG calculation. 

the ID phase in the present S~2 model. Furthermore, 
performing the level spectroscopy 10,11 ' (LS) analysis of 
the numerical results of exact-diagonalization calcula- 
tions, Nomura and Kitazawa 11 ' showed in the case of 
A = 1 that, with the increase of D from zero, the ground 
state changes from the Haldane state to the XYl state 
at D c i — 0.043 and after that from the XYl state to the 
LD state at D c2 — 2.39. Thus, it has been considered for 
a long time that the ID phase does not exist in the phase 
diagram of this model. 
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Fig. 3. (Color) The ground-state phase diagram on the A-D 
plane determined in the present work; in (b), part of (a) is en- 
larged. Note that the H and LD states belong to the same phase. 
See the text for more details. 



The purpose of the present work is to determine as pre- 
cisely as possible the ground-state phase diagram on the 
A-D plane of the 5 = 2 quantum spin chain described 
by the Hamiltonian (1). We mainly employ numerical 
methods based on the exact-diagonalization calculation 
for finite-spin systems with up to A = 12, the details of 
which are described below. For simplicity, we restrict our- 
selves to the A > and D > case. The full phase diagram 
will be discussed in our future reports. 

Figure 3 shows our final result for the ground-state 
phase diagram. In this phase diagram five states, that is, 
the XYl, LD, ID, Haldane and Neel states, appear as 
the ground states. As we will discuss later in detail, the 
Haldane and LD states belong to the same phase, which 
we denote as the Haldane/large-Z? (Haldane/LD) phase. 
Therefore, the phase diagram consists of four phases; 
they are the XYl, LD, Haldane/LD and Neel phases. 

In the following discussions for the determination 
of phase boundary lines, we denote, respectively, 
by E (N,M;pbc) and Ei(N, M ; pbc) the lowest and 
second-lowest energy eigenvalues of the Hamiltonian 
(1) with periodic boundary conditions within the sub- 
space determined by A and M. We also denote by 
Eo(N, M, P; tbc) the lowest energy eigenvalue of the 
Hamiltonian (1) with twisted boundary conditions, 
where an antiferromagnetic bond is twisted, within the 
subspace determined by A, M and P. Here, P(=+l or 
— 1) is the eigenvalue of the space inversion operator with 
respect to the twisted bond. 

Both the transition between the Haldane and ID 
states and that between the ID and LD states are 
the Gaussian transition, and those between the XYl 
state to one of the above three states are the 
Berezinskii-Kosterlitz-Thouless transition. 12 ' 13 ' Follow- 
ing the LS method, 10 ' 11 -' we should compare three 
excitation energies, Eq(N, 0, +1; tbc) — E (N, 0;pbc), 
E (N, 0,-1; tbc) - E (N, 0; pbc) and E (N, 2; pbc) - 
Eq(N, 0; pbc) in the A— i>oo limit. Namely, the ground 
state is one of the Haldane/LD, ID and XYl states, de- 
pending on whether the first, second or third excitation 
energy is the lowest among them. A physical and intu- 
itive explanation for this method is given in our recent 
paper. 14 ' 

From the above arguments, for a fixed D, the critical 
value of the Haldane- ID transition, Ac 11 ' 10 -*, and that of 
the ID-LD transition, Ac ID ' LD \ can be evaluated by the 
A^oo extrapolation of A<; H ' ID) (A) and A^ ID ' LD) (A), 
respectively, determined from the equation 10 ' 11 ' 

£ (A, 0,-1-1; tbc) = E (N, 0,-1; tbc) < E (N, 2; pbc). 

(2) 

Our calculations show that eq. (2) has two so- 
lutions A< H ' ID) (A) and A< ID ' LD) (A) for A = 8, 
10 and 12 in the region of 1.0<Z?<2.05, where 
A< H ' ID) (A)> A< ID ' LD) (A). In this region of D, we have 
estimated Ac H ' ID ^(A) and Ac ID,LD ' ) (A) by numerically 
solving eq. (2), and then extrapolated these results 
for A = 8, 10 and 12 to A— >oo by assuming that 
their A-dependences are quadratic functions of A~ 2 . 
Examples are shown in Fig. 4(a) and Table I. Thus, we 
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have obtained the Haldanc-ID and ID-LD boundaries, 
which are shown, respectively, by the red and blue solid 
lines in Fig. 3. 
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Fig. 4. (a) Examples of the Haldanc-ID and ID-LD transitions 
with £> = 1.82 and N=12, where A£j bc (Af) = Eo(N, 0, —1; tbc) 
-E (N,0, +l;tbc). From this, we obtain A<; H ' ID) (12) = 2.23793 
and Aj; H,ID) (12) = 2.16639 by the interpolation. We note that 
Eo(N, 2; pbc) is higher than Eo(N, 0, ±1; tbc) in this regime, (b) 
An example of the XY1-LD transition with D = 1.82 and N = 12, 
where AE£ _t (.rV) = E (N, 2; pbc)-E (N, 0, +1; tbc). From this, 
(XVI ID) 

we obtain A c ' (12) = 2.14265 by the interpolation. We 
note that E (N, 0, -1; tbc) is higher than E (N, 0, +1; tbc) in 
this regime. 



Table I. Examples of critical values of A in the case of D ■ 
obtained by the LS method. 
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We have found that the maximum value of D for 
which eq. (2) has two solutions exists for each value 
of TV and that, although this maximum value in- 
creases as TV increases, it approaches a finite value 
in the limit of N^oo. For example, eq. (2) has 
two solutions only when D<2.05, D < 2.14 and 
D < 2.17, respectively, for 7V = 8, 10 and 12. We have 
also found that, when eq. (2) has no solution, the 
relation E (N, M, +1; tbc) <E (N, M, -1; tbc) always 
holds. These show that, for sufficiently large values of 
D, the ID state does not appear in the ground-state 
phase diagram and the crossover between the Haldane 
and LD states may take place, which means that the 
Haldane and LD states belong to the same phase. In 
fact, the extrapolated lines of the phase boundary line 
between the Haldane and ID states and that between 
the ID and LD states, which are shown, respectively, by 
the red and blue dotted lines in Fig. 3, merge at the point 
(A, D) -(2.64, 2.19). 

The fact that the Haldanc and LD states belong to the 
same phase is very reasonable for the following reasons: 
(A) as discussed above, P= +1 for both states; (B) in the 
case in which open boundary conditions are imposed, it 
is clear that there exists no edge state in the LD state 
(see Fig. 1(b)), and also Pollmann et al. 15 ' 16 ) showed 



very recently that two 5=1/2 spins left at each edge 
in the Haldane state are in the (St t, S* ot ) = (1, 0) state, 
which leads to the no-edge state, at least if A > 1 and 
D>0 are satisfied; (C) Pollmann et al. 15 ) constructed a 
one-parameter matrix product state which interpolates 
the Haldane and LD states without any quantum phase 
transition. Although the Haldane and LD states are ap- 
parently different from each other from the valence-bond 
pictures in Figs. 1(a) and (c), they belong to the same 
phase. Here, we give two examples for such a situation. 
One is an antiferromagnetic 5=1 chain with the bond 
alternation and the on-site anisotropy, 17, 18 ' in which the 
LD state and the dimcr state belong to the same phase. 
The other is an antiferromagnetic 5=1 two-leg ladder, 19 ) 
where the Haldane state and the rung-dimer state belong 
to the same phase called the four-site plaquette singlet 
phase. 

In a similar way, for a fixed D, the critical values 
A C XY1 < H) andA (xyi < LD) 



of the ATFl-Haldane and XY1- 



LD transitions, respectively, can be obtained from 

E (N, 0,+l; tbc) = £ (^,2;pbc) < E (N, 0, -1; tbc). 

(3) 

An example is shown in Fig. 4(b) and Table I. For the 
critical value A c Jfyi ' ID ^ of the XYl-lD transition, we 
have to solve 

E (N, 0,-1; tbc) = E (N, 2; pbc) < E (N, 0, +1; tbc) 

(4) 

instead of eq. (3). Thus, we have obtained the phase 
boundary lines between the XY1 state and one of the 
Haldane, ID and LD states, which are depicted, respec- 
tively, by the brown, black and cyanic lines in Fig. 3. 

Finally, the phase transition between the Haldane and 
Neel states is of the 2D Ising type, since the symmetry 
is broken in the Neel state. In this case the phenomcno- 
logical renormalization group (PRC) method 20 ^ is useful 
for determining the phase boundary line. We have nu- 
merically solved the PRG equation, 

(/V -2){E 1 (N- 2,0; pbc) - E (N - 2,0; pbc)} 
= N{Ei(N, 0; pbc) - E (N, 0; pbc)} , 

to obtain the solution Ac H ' N ^(iV) for a given value of 
D. Then, we have extrapolated, to estimate the critical 
value A C H ' N ^ , these results for N = 8, 10 and 12 to N — »■ oo 
by assuming that their ./V-dcpendences are quadratic 
functions of (N — 1)~ 2 , and have determined the phase 
boundary line as shown by the green line in Fig. 3. It 
is expected that D c H ' N - ) approaches D^''^ — A in the 
A — > oo limit. This is because in the case of the Ising limit, 
governed by the Hamiltonian Hi s i ng — A ^-515^ 

\2 



(5) 



D ^2j(Sj) 2 , when D > A the ground state of the present 
system is the S| = state, and when D < A it is the Neel 
state with S| 



2(-iy and that with S? = 2(-1)J +1 state 
which are degenerate with each other. In the A — > oo 
limit the transition line between the Neel state and the 
Haldane state becomes of the first order, as is expected 
from the above Ising limit Hamiltonian Hising- Therefore, 
on the Haldane-Neel transition line, there should exist a 
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special point at which the transition changes from the 
second order to the first order as A increases. The de- 
tails of this point are left for a future study. 

The LS method 10, is based on the sine-Gordon the- 
ory or, equivalently, the c = 1 conformal field theory with 
the perturbation, where c is the central charge. Then, in 
order to check the applicability of the LS method, we 
have calculated the central charge c from the finite-size 
correction of the ground-state energy per one spin, 21-23 ) 

£ (A,0;pbc) £ o (A,0;pbc) _ ttv^c 

N ~ n^oo N 6N 2 { ' 

where v s is the spin wave velocity. The central charge c 
is defined only in the gapless region, and it becomes zero 
in the gapped region if the above formula is formally ap- 
plied. Along the line D = 1.9, where the ground-state 
changes as XY1 LD => ID Haldanc Necl with 
the increase of A from zero, we found that the values 
of "apparent c" are c~l for 0<A<2.7 and c~0 for 
A>2.7. This result shows that our system is effectively 
described by the sine-Gordon theory, which ensures the 
applicability of the LS method and is consistent with the 
smallncss of gaps in the LD, ID and Haldane states along 
this line. Nomura and Kitazawa 11 ) also found that c~l 
at the Haldanc-AYl and XY1-LD transition points on 
the A = 1 line and at the XY 1-Haldane transition point 
on the D = line. 

In conclusion, we have determined rather precisely the 
ground-state phase diagram on the A-D plane of the 
S = 2 quantum spin chain described by the Hamiltonian 
(1), employing various numerical methods based on the 
exact-diagonalization calculation. The resulting phase di- 
agram (see Fig. 3) consists of the XY1 phase, the Hal- 
dane/LD phase, the ID phase and the Neel phase. It is 
noted that in the Haldane/LD phase, the crossover be- 
tween the Haldane and ID states may occur. In our opin- 
ion, the most significant result of the present work is that 
the ID phase appears in the phase diagram. The existence 
of this phase was predicted by Oshikawa 6 ^ about twenty 
years ago. After then, considerable efforts by performing 
DMRG calculations 7-9 ^ were devoted to find this phase 
in the S = 2 case, but fruitful results were not obtained. 
We now think that the reason for this is clear. In the 
DMRG method, the distinctions between the XY1 state 
and the gapped states are made by extrapolating the ex- 
citation gaps of finite- N systems to N— >oo. Namely, the 
state is determined to be gapless or gapped depending 
on whether the extrapolated gap is zero or finite. Since, 
of course, the values of extrapolated gap have numer- 
ical errors, it is very difficult to discriminate between 
the gapless state and the state with a very small gap. 
This situation is clearly seen in Fig. 7 of the paper of 
Aschauer and Schollwock. 7 ) In other words, the appar- 
ent boundaries between the XY1 phase and the gapped 
phases determined by the DMRG method shift into the 
gapped region, because the system with £ > N behaves to 
be gapless, where £ is the correlation length. In fact, in 
the case of A = l, the transition point between the XY1 
phase and the LD phase is estimated as _D c2 ^3.0 by 
the DMRG method, 9 ) which is considerably larger than 



D c2 = 2.39, n ) obtained by the LS method. As can be 
seen from Fig. 3, the region of the ID phase is fairly nar- 
row, and thus the whole region of the ID phase was com- 
pletely masked by the "apparent XY 1 region" . Although 
we have tried to detect the edge state with S z = 1/2 (see 
Fig. 1(b)) around each end of the open chain in the ID 
state by the DMRG method, this has also led to fruitless 
results. The reason why the LS method 10 ' n ) is very use- 
ful is undoubtedly the fact that it can sharply distinguish 
three phases, the XY1, ID and Haldane/LD phases, by 
the crossing of the excitations under periodic and twisted 
boundary conditions. 
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